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Abstract. Interaction properties of complex solitons are studied for the two J7(l)-invariant 
integrable generalizations of the mKdV equation, given by the Hirota equation and the Sasa- 
Satsuma equation, which share the same travelling wave (single-soliton) solution having 
a sech profile characterized by a constant speed and a constant phase angle. For both 
equations, nonlinear interactions where a fast soliton collides with a slow soliton are shown 
to be described by 2-soliton solutions that can have three different types of interaction 
profiles depending on the speed ratio and the relative phase angle of the individual solitons. 
In all cases the shapes and speeds of the solitons are found to be preserved apart from a 
shift in position such that their center of momentum moves at a constant speed. Moreover, 
for the Hirota equation, the phase angles of the fast and slow solitons are found to remain 
unchanged, while for the Sasa-Satsuma equation, the phase angles are shown to undergo a 
shift such that the relative phase between the fast and slow solitons changes sign. 



1. Introduction 

Solitons are solitary waves that retain their shape and speed after undergoing collisions or 
other nonlinear interactions [1]. A widely familiar example is the sech!^ solitary wave solution 
of the Korteweg-de Vries (KdV) equation ut + 6uUx + Uxxx = 0. This solution u{t,x) = 
(c/2)sech2 {^{x - ct)/2) is a stable travelling wave that is single-peaked and uni-directional. 
It carries mass, momentum, energy, as well as Galilean momentum (associated with the 
motion of center of mass), which are constants of motion for the KdV equation. A collision 
occurs when a faster solitary wave overtakes a slower solitary wave. (See the animations 
at http://lie.math.brocku.ca/~sanco/solitons/kdv_solitons.php ) Remarkably, the only net 
effect of the collision is that the faster wave is shifted forward in position while the slower 
wave is shifted backward in position, where these shifts depend solely on the speeds of both 
waves and do not affect the center of mass of the two waves (which moves at constant speed 
throughout the collision). 

During a collision, KdV solitary waves interact nonlinearly such that [2, 3] their peaks 
either first merge together and then split apart if the speed ratio of the waves is greater than 
3 or first bounce and then exchange shapes and speeds if the speed ratio of the waves is less 
than 3. There are many different alternative ways to interpret this nonlinear interaction, 
such as the faster wave getting stretched while the slower wave is squeezed underneath it, 
or as the faster wave emitting an intermediate wave that is absorbed by slower wave. (See 
Ref. [4] for a comprehensive survey of interpretations.) The same interaction properties also 
hold more generally for pair-wise collisions of any number of KdV solitary waves. 
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Solutions of the KdV equation describing collisions of n > 1 solitary waves are called 
n-solitons and have been obtained by many different methods (e.g. auto-Backlund trans- 
formations, nonlinear superposition formula, Hirota bilinear equations, inverse scattering, 
dressing equations), all of which rely on the underlying integrability of the KdV equation, 
in particular the existence of a Lax pair [2]. An important consequence of this integrability 
is that KdV solitons have constants of motion consisting of mass, momentum, Galilean mo- 
mentum, and energy, which are defined by conserved integrals [5] involving t,x,u,Ux, plus 
an infinite number of higher-order "energies" (involving higher order x-derivatives of u) [1] . 

In this paper, we study the coUision properties of sohtons of the modified KdV (mKdV) 
equation 

Ut + au^u^ + /3uxx^ ^ (1.1) 
and its integrable C/(l)-invariant generalizations 

ut + a\u\'^Ux + Puxxx ^ (1.2) 
ut + la{uu^ + 3u^u)u + I3u^^^ ^ (1.3) 

where u denotes the complex conjugate of u, and |m| denotes the modulus of u. Here a 
and (5 are arbitrary positive constants. These two generalizations are known to be [6, 7] the 
only complex versions of the mKdV equation that possess a Lax pair with the same scaling 
symmetry 

t XH, X Xx, X~^u (1.4) 

admitted by the mKdV equation (1.1). Both generalizations also have an additional U{1) 
phase symmetry 

u — )■ exp{i(f))u (1.5) 

In this form, all three equations (1.1), (1-2), (1.3) share the same solitary wave solution (i.e. 
a 1-soliton) 




u{t,x) — ±\ — sech. [ . —(x — ct) ] (1.6) 



where c > is the wave speed. As shown by the results in Ref. [8] on constants of motion for 
equations of complex mKdV form, the solitary wave solution (1.6) has mass, momentum and 
energy, which are given by counterparts of the KdV conserved integrals involving just UjU^] 
in addition, although there is no counterpart of KdV Galilean momentum for this solution, 
it has an analogous Galilean energy given by a conserved integral which is related to the 
motion of center of momentum. 

In the case of the real mKdV equation (1.1), the 1-soliton (1.6) has an up or down orien- 
tation corresponding to the plus or minus sign of u, which comes from the discrete reflection 
symmetry u —u of this equation. Thus, there are two different types of real mKdV soliton 
colhsions, where (up to reflection) the fast and slow sohtons in the coUision have either the 
same orientation or opposite orientations. 

In both cases of the complex mKdV equations (1.2) and (1.3), the sign of the 1-soliton 
(1.6) can be absorbed into an arbitrary constant phase 

/6c / ic \ 

— exp(i0)sech I . / — (x — ct) ] , (j) — const. (1.7) 
a \y P J 
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due to the U{1) symmetry (1.5). Consequently, collisions of two complex mKdV solitons 
(1.7) involve a relative phase angle, given by the difference of the phase angles of the fast 
and slow solitons in the collision. This relative phase therefore parameterizes the types of 
collisions. An interesting question we will study is whether the two phase angles are altered 
in a collision, i.e. do the fast and slow solitons each get shifted in phase as well as in position? 

In Sec. 2, the 2-soliton solution that describes collisions of fast and slow solitary waves for 
the mKdV equation (1.1) is reviewed. The properties of these collisions depend only on the 
ratio of speeds and the relative orientation of the two waves. In particular, we show that if 
the waves have the same orientation then their nonlinear interaction in a collision either is a 
merge-split type when their speed ratio is greater than the value (7 + 3\/5)/2, or otherwise is 
a bounce-exchange type when their speed ratio is less than this critical value (7 + 3v^)/2. In 
contrast, if the waves have opposite orientations then their nonlinear interaction instead is a 
completely different type in which (regardless of their speed ratio) the slow soliton gradually 
is first absorbed and then emitted by the fast soliton. In all three interactions, we show that 
the net effect of the collision is solely that the faster wave is shifted forward in position while 
the slower wave is shifted backward in position, such that the center of momentum of the 
two waves is unaffected (moving at a constant speed throughout the collision). 

In Sec. 3, we consider the complex mKdV equation (1.2), which is commonly called the 
Hirota-mKdV equation [9]. We carry out an asymptotic analysis of the 2-soliton solution 
describing collisions of fast and slow solitary waves. Our results show that the interaction 
properties of the two waves depend only on their speed ratio and their relative phase angle 
such that a bounce-exchange type of interaction occurs when the speed ratio is less than 
a critical value given by a certain explicit function of the relative phase angle and that 
otherwise a merge-split or absorb-emit type of interaction occurs depending on whether the 
relative phase angle is less than or greater than a certain critical value in terms of the speed 
ratio of the waves. For each type of interaction, we find that the phase angles of the two 
waves remain unchanged and the only effect of the collision is to produce a respective forward 
and backward shift in the positions of the faster and slower waves. In particular, these shifts 
are found to depend only on the speeds of the two waves, but not on their phases angles, 
such that the center of momentum of the two waves moves at a constant speed throughout 
the collision. 

In Sec. 4, we consider the other complex mKdV equation (1.3), which is known as the 
Sasa-Satsuma-mKdV equation [10]. We write down the 2-soliton solution in an explicit form 
parameterized by the speeds and phase angles of the fast and slow solitary waves (which has 
not appeared previously in the literature). Through an asymptotic analysis of this solution, 
we find the interesting new result that the phase angles of the two waves in a collision undergo 
a shift such that the relative phase angle changes sign. In addition, the positions of the fast 
and slow waves display a respective forward and backward shift which depends on both the 
speeds and the relative phase angle of the waves. We show that these position shifts preserve 
the center of momentum of the two waves in the collision, while the phase- angle shifts are 
related to an invariance property of 2-soliton solution with respect to space-time refiection 
combined with phase conjugation. Finally, we also derive the detailed interaction properties 
of the two waves. We show that the waves exhibit a bounce-exchange type of interaction 
only when their the relative phase angle is less than arccos(— 1/6) ~ O.SStt and their speed 
ratio is less than a critical value given by a certain explicit function of the relative phase 
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angle (which is different than the function arising for the Hirota-mKdV 2-soliton sohition). 
For any speed ratio greater than this value, or for any relative phase angle greater than 
arccos(— 1/6) ~ O.SStt, we find that the waves exhibit a merge-split type of interaction if 
their relative phase angle is less than a certain critical value in terms of the their speed 
ratio (which is again different than the critical angle found for the Hirota-mKdV 2-soliton 
solution), and that otherwise the waves exhibit an absorb-emit type of interaction. 

Last, some features of the soliton collisions for the Hirota and Sasa-Satsuma equations are 
compared in section Sec. 5. 

In the appendix, we provide a short derivation of the 2-soliton solution for the Sasa- 
Satsuma-mKdV equation. Hereafter, by scaling variables, we will put 

a = 24, /3 = 1 (1.8) 

for convenience. 



2. REAL MKdV soliton COLLISIONS 

For the mKdV equation 

ut + 24:U^u^ + u^^^ = (2.1) 

we first recall the conserved integrals defining counterparts of KdV mass, momentum, and 
energy. These integrals are given by [11] 

f +00 

M 



/+00 
u dx (2.2) 
■oo 

dx (2.3) 



oo 

+ 0O 



oo 
hoo 



j ^ul-2uUx (2.4) 

Although there is no counterpart of KdV Galilean momentum, the mKdV equation has the 
extra conserved integral [11] 

C= t{^ul - 2u') - x-v? dx (2.5) 



which defines a Galilean energy related to center of momentum as given by 



1 /■"^°° E 
X{t) = ^ / XV? dx^ A'(O) + 6-i (2.6) 

' J —oo 



where 

C^tS- -VX(t) = C(0) = --VX(0) (2.7) 
6 6 

These integrals (2.2)-(2.5) are constants of motion for all smooth solutions u{t,x) with 
asymptotic decay it — > as x — > ±oo. In particular, this includes the solutions describing 
solitary waves and their collisions. 

The 1-soliton for the mKdV equation is given by [12] 

u{t, X) = ±-sech (V^O = ^ + exp(2VH0 ^'"'^ 
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with speed c > and up/down orientation s = ±1, where 



i^x-ct (2.9) 

is a moving coordinate. This solution describes a stable travelling wave that is single-peaked 
and uni-directional. Its height relative to = is ±-^/2, and its width is proportional to 
1 / ^/c. It has constants of motion 

M^s\, V^^, S^^, C = 0, (2.10) 
while its center of momentum is 

Xit) = 6^t = c± (2.11) 

which coincides with the position of the peak. The initial position of the wave can be shifted 
arbitrarily by means of a space translation x ^ x — xq applied to the moving coordinate 
(2.9), so then ^ — x — ct — xq and X{t) = xq + ct. This changes the Gahlean energy of the 
wave to be 

C = -^0^0/12, (2.12) 
while the mass, momentum and energy are unchanged. 

2.1. 2-soliton solution. 

Collisions where a fast soliton with speed Ci overtakes a slow soliton with speed 
C2 are described by the well-known 2-soliton solution [12] which depends on the ori- 
entations Si and 6"2 of the respective solitons. (See the animations of collisions at 
http://lie.math.brocku.ca/~sanco/solitons/mkdv_solitons.php ) We will write this solution 
in the rational-exponential form 

u{t,x)^j (2.13) 

given by 

G =K(siv/^cxp(v/^ei)(l + exp(2v/S^6)) + S2v^cxp(v/S^6)(l + exp{2^^i))) (2.14) 
F=l + 2siS2(k' - 1) cxp(v/Siei + ^M^2) + «;'(exp(2V^ei) + exp{2^^2)) 

+ exp(2(V^ei + 7^6)) (2.15) 



with 



where 



Sl = ±l, S2 = ±l, K= > 1 (2.16) 



^1= X - Cit - Xi, = X - C2t- X2 (2-17) 



are moving coordinates centered at initial positions x = Xi and x = X2 respectively. 

We now examine the asymptotic form of the 2-soliton solution (2.13)-(2.15) as t — )■ ±oo. 

To proceed, we consider = const, and express ^2 — Ci+C terms of ^1 and C = tAc—Ax, 
where Ac = ci — C2 > is the relative speed of the moving coordinates, and Ax — X2 — xi 
is the separation of the centers of the moving coordinates. Note t ±00 corresponds to 
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C — )■ ±cxo. We then asymptotically expand F and G for large C, with held fixed. This 
yields, after neglecting subdominant exponential terms, 

Q ^ Hi«^v^exp((^ + 2^)^i) exp(2^C) C +00 18) 
\siK^exp(^^i) C,^-Qo 

F r^l^^^ exp(2y^,^i) + exp(2(^ + v^)6)) exp(2y^C) C ^ +00 -^gx 
~ |l + K2exp(2^ei) 

and hence 



^ {\fci/i^) cxp(ycrci) ^ ^ 

"-{"'^ilCfir/^if^^'* (2.20) 



y/clKexp(-v/cI^ij 
'^l + «:2exp(2y5r6) °° 

Thus, in this expansion the 2-soliton solution asymptotically reduces to the form of a 1- 
soliton solution 

u{t, x) ~ gi f^^'^'^lf^l}^ = uf, t^ioo, 6 = const. (2.21) 
1 + exp(2^^i ) 

in terms of a moving coordinate 

^1 =^i-af (2.22) 

which is shifted relative to ^1 by 

af = ±ln(K)/^ (2.23) 

To continue, we next consider ^2 = const, and express ^1 = S ~ C iii terms of ^2 and 
( = tAc — Ax again. By asymptotically expanding F and G for large C with ^2 held fixed, 
and neglecting subdominant exponential terms, we obtain 

^ _ |s2/tv^exp(yS^G) ^2 24) 

\s2/«V^exp((y^+ 2v^)^2) exp(-2y^C) ( -00 

^ ^ f 1 + «2 exp(2^e2) C ^ +00 

1 exp(2^e2) + exp(2(^ + ^51)6)) exp(-2^C) C ^ -00 



and thus 



/o7/vexp(yr7(2) ^ ^ , 

C +00 



u 



■l + /^^exp(2^6) (2.26) 
(yc^/ft^) exp(^^2) A >_ 

'l + (l//^2)exp(2^e2) ^ 
This asymptotic expansion of the 2-soliton solution again has the form of a 1-soliton solution 



C2exp(^^2 



u{t, x) ~ ^ =U2, ±00, C2 = const. (2.27) 

1 + exp(2^^f ) 

in terms of a moving coordinate 

^2"= = 6 - «^ (2.28) 
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which is shifted relative to ^2 by 

af = =Fln(K)/v/c^ 



(2.29) 



2.2. Constants of motion and asymptotic position shifts. 

The expansions (2.21) and (2.27) show that for t — > ±00 the 2-soliton solution (2.13)- 
(2.15) asymptotically has the form of a superposition u ^ + uf oi a fast soliton uf with 
speed ci and up/down orientation si and a slow soliton uf with speed C2 and up/down 
orientation S2, whose positions are determined by the moving coordinates and ^f. (This 
result generalizes the well-known analysis in Ref. [13] which considered the case si = S2 = 1.) 
See Fig. 1 and Fig. 2. 
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Figure 1. mKdV equation 2-soliton interaction with Ci/c2 = 2, Si/s2 = 1 
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Figure 2. mKdV equation 2-soliton interaction with Ci/c2 = 2, Si/s2 



In the asymptotic past (t — )■ —00), = gives the position of the peak of the fast 
soliton, while in the asymptotic future (t +00), the position of the peak is instead given 
by = 0. These asymptotic positions lie on straight lines in space-time 

x = xi + af + cit (2.30) 
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Comparing the asymptotic past with the asymptotic future, we see that the fast sohton 
retains its shape and speed but gets shifted forward in position as given by 

Ax, = a+ - ar = ^In f ^+ ^ ] > Q (2.31) 

Similarly, the slow soliton retains its shape and speed, while the position of its peak in 
the asymptotic past and future lies on the straight lines x = X2 + 02 + C2t. So we see that 
the slow soliton gets shifted backward in position as given by 

Ax2 ^at-a2^-^ln( ^^±v^^ < (2.32) 



The asymptotic shifts (2.31) and (2.32) do not depend on the orientations of the fast and 
slow solitons. Moreover, these shifts satisfy the relation 

^Axi + y^Ax2 = (2.33) 

which can be understood as a consequence of the motion of the center of momentum of the 
2-sohton solution (similarly to the same result known for the KdV 2-soliton solution [14]). 

In particular, because ~ + is a superposition as i — >■ ±00, the conserved mass, 
momentum and energy of u are given by 

M = A4i + M2 = (si + S2)^ (2.34) 



V^Vi + V2= \ ^""^ (2.35) 

S = S, + 82= + (2.36) 

in terms of the mass, momentum and energy individually associated with the fast and slow 
solitons and u^- Similarly, the conserved Galilean energy of u is given by 

(^^i + %) + (^^2 + 02 ) _ xi + at)^i + {x2 + a^)^ , . 

C-C1+C2- (2.37) 

which simplifies to 

^_ ^l^ + ^2v/^_^^_lp;t.(^) (2.38) 
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due to the relation af ^/ci + a^y^ = 0. Thus the center of momentum of u is 

^ V^C'-i + ^'iO+ v^C'-2 + C2t) ^ PiAf(f) + P2^y^(0 .2 39) 

where Xi{t) — x, -\- af -\- Cit and X^{t) — X2 + af + C2t are the respective centers of 
momentum of the fast and slow solitons in the asymptotic past and future. Therefore, we 
see that the center of momentum of the 2-soliton solution moves at constant speed 

and consequently the asymptotic shifts Axi and Ax2 in the positions of the two solitons 
and uf are constrained to satisfy ViAxi + V2AX2 — which explains the relation (2.33). 



2.3. Interaction profile. 

We now study the interaction profile of the 2-soliton solution (2.13)-(2.15), which we will 
write in the equivalent form 

_ K{siy/c^ cosh{92) + g2A/c2 cosh(6'i)) 

' ~ siS2{k^ -1) + k^ cosh(^i - 62) + cosh(^i + 62) ^' ' 

in terms of 

Oi = v^6, O2 = (2.42) 

For simplicity, by means of suitable time and space translations t ^ t — to, x ^ x — xq, we 
shift the centers of the moving coordinates (2.17) to the positions 

xi = xo - cito = 0, X2 = xo - C2to = (2.43) 

Then the resulting 2-soliton solution (2.41) is invariant under a combined space-time reflec- 
tion X —X, t —t, and its center of momentum is simply X{t) = ct in terms of the speed 
(2.40). In the asymptotic past and future {t — )■ =Foo), this solution describes a superposition 
u uf + uf oi fast and slow solitons uf and whose centers of momentum are given by 
^I'i't) — af + cit and X^{t) — + C2t where 

AX^{t) = X^{t) - X^{t) = ±^{Axi - Ax2) - (ci - C2)t (2.44) 

is the separation between the peaks of and uf in terms of the asymptotic shifts (2.31) 
and (2.32) for t ^oo. 

Because of the space-time reflection invariance of u, the separation between the fast and 
slow peaks in u will be a minimum at time t = when u is an even function of x. Qualitatively 
speaking, t — will be the moment of greatest nonlinear interaction between the fast and 
slow solitons. We will therefore refer to the shape of w at i = as the interaction profile of 
the 2-soliton solution. Since this profile u{0, x) is even in x, its shape can be characterized 
by the convexity 

W(0,a;)a;x|^^0 = ^(Ci - C2){^/C^- \/c^){y/c^- {Siy/c^ - S2y/c^f) (2.45) 

Note that the sign of the convexity (2.45) depends only on the ratio of speeds and the relative 
orientation 

r = ci/c2, s = si/s2 (2.46) 
of the fast and slow solitons. In particular, since ci > C2 > 0, we have 

a = sgn(M(0,x)^4^J = sgn((2s + l)V^-r-l) (2.47) 

with r > 1 and s = ±1. 

The case where the two solitons have the same orientation is given by the relation s = 1. 
In this case, the convexity sign is determined by 



a = sgn(3Vr - r - 1) = sgn(|(3 + -\/5) - ^/r) 
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(2.48) 



via factorization of (2.47) and use of the inequality r > 1. Hence the sign is indefinite such 
that 

{<0 if r > (7 + 3v/5)/2 
>0 if r < (7 + 3y5)/2 (2.49) 
= if r = (7 + 3v^)/2 

This result implies that the interaction profile ^(0, x) will have either a single peak at x = 
if ci/c2 > (7 + 3v^)/2, or a double peak around x = if Ci/c2 < (7 + 3v^)/2. These peaks 
will be positive or negative depending on the sign of si (= S2) and will have an exponentially 
diminishing tail. In the case of a single peak, the fast and slow solitons interact by first 
merging together at t = x = and then splitting apart, while in the case of a double peak, 
the fast and slow solitons interact by first bouncing and then exchanging shapes and speeds 
at t = a; = 0. We will call these cases, respectively, a merge-split and bounce- exchange 
interaction. See Fig. 3 and Fig. 4. 
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Figure 3. mKdV equation 2-soliton interaction with Ci/c2 = 9, Si/s2 = 1 
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Figure 4. mKdV equation 2-soliton interaction with Ci/c2 = 4, Si/s2 = 1 
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The value Ci/c2 = (7 + 3-\/5)/2 which separates these two types of interaction profiles will 
be called the critical speed ratio. 

The other case, where the two solitons have opposite orientations, is given by s = — 1. In 
this case, the convexity sign (2.47) is strictly negative for all r > 1, 

cr = -sgn(v/r + r + 1) < (2.50) 

implying that the interaction profile u{0,x) will have a positive or negative peak at x = 
depending on the sign of si (= — S2). Since S2 has the opposite sign, the profile will also have 
a pair of negative or positive side peaks around a; = 0, with an exponentially diminishing tail. 
The interaction between the fast and slow solitons in this case consists of the slow soliton 
gradually being first absorbed by the front side of the fast soliton and then emitted from the 
back side of the fast soliton. We will call this an absorb-emit interaction. See Fig. 5. 
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Figure 5. mKdV equation 2-soliton interaction with Ci/c2 = 9, Si/s2 = — 1 



3. Hirota-mKdV soliton collisions 
For the Hirota equation 

+ 24|m|\^ + M^^^ = (3.1) 

we note, firstly, there is no conserved integral for mass since neither J^°° u dx nor \u\ dx 
is a constant of motion. Secondly, we recall the conserved integrals for momentum, energy, 
and Galilean energy are given by [8] 

/+00 
\u\^ dx (3.2) 
-00 

r+oo 

£= |y^|2_2|u|^cia; (3.3) 



-00 
+00 



2' 



1 1 

C= I t{-\u^\^ -2\u\'^) - x-\u\^ dx (3.4) 
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which are constants of motion for all smooth solutions u{t,x) with asymptotic decay u ^0 
as X ^ ±00. These integrals are related to the center of momentum 

1 /"+°° £ 

X{t) = - x\u\'^ dx = A'(O) + 6-i (3.5) 
' J —00 I 

since 

C = tS- lvX(t) = C(0) = -lvX{0) (3.6) 
6 6 

This is the same relation that holds for the mKdV constants of motion. 
The 1-soliton solution for the Hirota equation is given by 

with speed c > and phase — tt < < tt, where 

^ = x-ct (3.8) 

is a moving coordinate. This solution describes a stable uni-directional travelling wave whose 
amplitude |m| is the same as the amplitude of the mKdV solitary wave (2.8). Therefore, its 
constants of motion (3.2)-(3.4) and its center of momentum (3.5) are also the same as those 
for the mKdV solitary wave. In particular, the position of the peak amplitude coincides with 
the center of momentum X{t) = ct, which can be shifted arbitrarily by means of a space 
translation 

X ^ X — xq (3.9) 

applied to the moving coordinate (3.8). 

3.1. 2-soliton solution. 

We now write down the 2-soliton solution of the Hirota equation describing collisions where 
a fast soliton with speed Ci and phase 0i overtakes a slow soliton with speed C2 and phase 02. 
(See the animations of collisions at http://lie.math.brocku.ca/~sanco/solitons/hirota.php ) 
This solution has the rational-exponential form (2.13) given by [9] 

G =K{y/c^exp{i(f)i + Vci{i)(l + exp(2^^2)) + y/c^exp{t(})2 + x/ciC2)(l + exp(2^{i))) 

(3.10) 

F =1 + 2cos(0i - 02)(/t' - 1) exp(v/3r6 + v^6) + n\exp{2^ii) + exp{2^^2)) 

+ exp(2(y^ei + v^6)) (3.11) 



in terms of 



where 



. = ^^±^>1 (3.12) 



= X — Cit — Xi, ^2 = X — C2t — X2 (3.13) 

are moving coordinates centered at initial positions x = Xi and x = X2 respectively. 

The asymptotic form of the 2-soliton solution u{t,x) = G/F ajst ^ ±00 is easily derived by 
the same moving-coordinate expansions considered for the mKdV 2-soliton solution. These 
expansions yield, after neglecting subdominant exponential terms, 

= «f , t^ico, 6= const. (3.14) 

l + exp(2^(^i -af)) 
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and 

u{t, x) ~ — — — = M2 , ±00, 6 = const. (3.15) 

1 + exp(2^(^2 - ^2 )) 

where af , are shifts in the moving coordinates ^1, ^2 given by (2.23) and (2.29) respectively. 

Thus, the 2-sohton solution (2.13) and (3.10)-(3.11) asymptotically has the form of a 
superposition -u ~ -uf + -u^ for t — )■ ±00, where is a fast soliton (3.7) with speed Ci and 
phase 01 and where is a slow soliton (3.7) with speed C2 and phase 02, whose positions are 
determined by the shifted moving coordinates (2.22) and (2.28). This solution has the same 
conserved momentum (2.35), energy (2.36), Galilean energy (2.38), and center of momentum 
(2.39) as the mKdV 2-soliton solution. 

3.2. Asymptotic position shifts. 

For the 2-soliton solution u of the Hirota equation, the asymptotic expansions (3.14) and 
(3.15) for t — )■ —00 compared to t — )■ +00 show that the fast soliton ui retains its shape, 
speed and phase, but gets shifted forward in position by 



Axi = a+ - ar = ^In ( v^±vg ) > (3.16) 

while the slow soliton U'z similarly retains its shape, speed and phase, but gets shifted back- 
ward in position by 

Ax2 = a+ - a2- = In f ^+^ ) < (3.17) 



These expressions are the same asymptotic shifts seen for the collision of mKdV solitons. In 
particular, the shifts (3.16) and (3.17) here are independent of the phases of the fast and 
slow solitons and also satisfy the center of momentum relations (2.33) and (2.44). 

3.3. Interaction profile. 

In the same way as for the mKdV 2-soliton solution, we will now write the Hirota 2-soliton 
solution given by (2.13) and (3.10)-(3.11) in the equivalent rational-cosh form 

_ ft:(^exp(i0i) cosh(6'2) + Vc^exp(i02) cosh(6'i)) 18^1 
^^'^>~ («;2 _ 1) cos(0i - 02) + C0Sh(^i - 62) + cosh(^i + ^2) 

with 

Oi = v^Ci, e.2 = (3.19) 
By shifting the centers of the moving coordinates (3.13) to the positions 

xi = xo - Clio = 0, X2 = Xq- C2to = 0, (3.20) 

we then see that the resulting solution (3.18) of the Hirota equation is invariant under the 
combined space-time reflection x — )> —x, t —t. As a consequence, this solution u describes 
a collision of fast and slow solitons and whose separation will be a minimum at time 
t — Q when the amplitude \u\ is an even function of x. This can be understood to be the 
moment of greatest nonlinear interaction between the fast and slow solitons in the collision. 
The shape of |m| at t = therefore defines the interaction profile of the 2-soliton solution 
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u. Since this profile is even in x, its shape is characterized by the convexity of |w(0, x)| at 
a; = 0. An explicit calculation of the convexity yields 

I /p, x| I _ (ycl - V^)^(4ciC2 -cl-4 + v/cIc^(Cl + C2) COS(01 - 02)) 

\u[y), - + e2 + 2V5I^cos(0i - ^2))^/^ ^ ^ ^ 

The sign of the convexity (3.21) depends only on the ratio of speeds and the relative phase 

r = Ci/c2, A(j) = 01 — 02 mod 27r (3.22) 

of the fast and slow solitons. In particular, since Ci > C2 > 0, we have 

a = sgn(|ii(0,x)U4^p) = sgn(4r - - 1 + V^(r + 1) cos A0) (3.23) 

with r > 1 and — tt < A0 < tt. This expression (3.23) is a quartic polynomial in ^/r. By 
factorizing 

4r - - 1 + Vr(r + 1) cos A0 = 

(|(cosA0+ a/cos^ A0 + 24) Vr - r - l)(|(cosA0+ -y^cos^ A0 + 24) Vr + r + l) 

and using the inequality 



a/cos^ A0 + 24 > ± cos A0 

we obtain 

a = sgn(|(cos A0 + ^cos^ A0 + 24) Vr - r - 1) (3.24) 
which is a quadratic polynomial in ^/r with two real positive roots 



ri= i(cosA0+ Vcos2A0 + 24) ± A/(i(cosA0+ Vcos2A0 + 24))2 - 1 >0 (3.25) 



These roots satisfy r+r_ = 1, where the case r+ = r_ = 1 occurs iff |A0| = tt. Hence we 
have 

r+ > 1 > r_ (3.26) 
in all cases. Thus the convexity sign (3.24) is determined by 

a = sgn( - v^) (3.27) 



which is indefinite such that 




(3.28) 



where 



r+ = 2 + i ( cos A0 + a/cos2 A0 + 24 ) ( cos A0 + y |(cos A0 + a/cos^ A0 + 24) cos A0 + 2 ) 

(3.29) 

The result (3.28) implies that the 2-soliton interaction profile |M(0,a:)| will have either a 
single peak at a; = if Ci/c2 > r_|_ or a double peak around a; = if Ci/c2 < r+, as determined 
by the critical speed ratio (3.29) in terms of the relative phase angle A0. 

In the case of a single peak, the profile has one of two different shapes depending on 
whether |A0| is greater than or less than a certain critical value given by some function 
of Ci/c2 that is determined by the conditions |m(0, a;)!-,;-,; = |M(0,a;)|a; = for existence of a 
saddle point at some x ^ (which we can solve for numerically). For |A0| below the critical 
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value, the shape of |m(0,x)| is simply a single peak with an exponentially diminishing tail. 
In this case the fast and slow solitons undergo a merge-split interaction, i.e. where they first 
merge together ai t = x = and then split apart: See Fig. 6 and Fig. 7. For |A0| above the 
critical value, the shape of |m(0,x)| consists of a pair of side peaks around the main peak at 
X = 0. In this case the fast and slow solitons undergo an absorb-emit interaction, i.e. where 
the slow soliton gradually is first absorbed by the front side of the fast soliton and is then 
emitted from the back side of the fast soliton. See Fig. 8 and Fig. 9. 




Figure 6. Hirota equation 2-soliton interaction with Ci/c2 = 3, — 02 = O.SSSvr 



N |u| |u| 




Figure 7. Hirota equation 2-soliton interaction with Ci/c2 = 50, 0i — 02 = O.Stt 
The interaction in the special case when |A0| equals the critical value is shown in Fig. 10. 

In contrast, in the case of a double peak, the profile |M(0,a;)| always has an exponentially 
diminishing tail, regardless of the relative phase angle |A0|. This case describes the fast and 
slow solitons undergoing a bounce-exchange interaction, i.e. where they first bounce and 
then exchange shapes and speeds at t = x = 0: See Fig. 11. 

The range of values of Ci/c2 and |A0| that characterize these three different types of 
interaction is shown in Fig. 12. 
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Figure 8. Hirota equation 2-soliton interaction with Ci/c2 = 3, 0i — 02 = O.SSSvr 




Figure 10. Hirota equation 2-soliton interaction with Ci/c2 = 3, 



>1 — V2 



0.855987r 
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Figure 11. Hirota equation 2-soliton interaction with Ci/c2 = 3, 0i — 02 = O.Stt 
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Figure 12. Types of 2-soliton interactions for the Hirota equation 
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4. Sasa-Satsuma-mKdV soliton collisions 



For the Sasa-Satsuma equation 

Ut + 6{uUx + 3Uj:U)u + Uxxx = (4.1) 

we begin by remarking [8] that it has the same conserved integrals as the Hirota equation for 
momentum (3.2), energy (3.3), and Gahlean energy (3.4), which define constants of motion 
for all smooth solutions u(t,x) with asymptotic decay m — )■ as x — )■ ±oo. In addition these 
integrals have the same relation to the center of momentum (3.5)-(3.6) that holds for the 
Hirota equation. 

The 1-soliton solution for the Sasa-Satsuma equation has the same form (3.7)-(3.8) as 

the Hirota 1-soliton with speed c > and phase — vr < (j) < n. We recall that this solution 
describes a stable uni-directional travelling wave u whose amplitude \u\ is the same as the 
amplitude of the mKdV solitary wave (2.8), and hence u also has the same constants of 
motion (3.2)-(3.4) and center of momentum (3.5) as those of the mKdV solitary wave. Thus, 
the position of the peak amplitude of u coincides with the center of momentum X{t) — ct 
which can be shifted arbitrarily via a space translation x ^ x — xq applied to the moving 
coordinate (3.8), yielding 

„(,.)^fexp,,,sech(V3.) = 4^^g*±f) ,4.2, 

with 

^ — X — ct — xq (4.3) 

4.1. 2-soliton solution. 

The 2-sohton solution of the Sasa-Satsuma equation describing collisions where a fast 
soliton with speed Ci and phase 0i overtakes a slow soliton with speed C2 and phase 02 has 
not appeared previously in an explicit form [15]. We give a simple derivation of this solution 
in the appendix, based on using a rational-exponential ansatz similar to the form of the 
2-soliton solution of the Hirota equation. This derivation yields 

u{t,x)^j (4.4) 

with 

G (^^/c^exp{i(f)l + ^/c^Ci) ((^/cl + ^/c^exp(2^(02 - 0i))) exp(2^^2) + 7) 

+ ^/c^exp{i4>2 + VC2C2)(7exp(2v^^i) + ^/cT exp(2i(0i - ^2)))^ (4.5) 

F =1 + 2 cos(0i - h){{V^ + v^)p' - 1) eMVcI^i + v^6) 

+ 7p2(exp(2v/^ei) + exp(2^6)) + exp(2(V^ei + v^6)) (4-6) 

in terms of 



p= -^f^f >0 (4.7) 

Vci- ^/C2 

and 

7 = |vciexp(i(02 - 0i)) + V^exp(i(0i - ^2))! = \Jci + C2 + 2^ciC2 cos 2(0i - ^2) > 

(4.8) 
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where 

^l^X — Cit — Xi, ^2^X — C2t — X2 (4.9) 

are moving coordinates centered at initial positions x = Xi and x = X2 respectively. (See the 
animations of collisions at http://lie.math.brocku.ca/~sanco/solitons/sasa-satsuma.php ) 

We will now examine the asymptotic form of the 2-soliton solution (4.4)-(4.6) as t ^ ±00 
by means of the same moving-coordinate expansions used for the mKdV 2-soliton solution. 

First we hold ^1 fixed and asymptotically expand F and G for large C = tAc — Ax with 
C2 — ^1 + C- This yields, after neglecting subdominant exponential terms, 

{-^ycrexp(i0i + (^+ 2^)^i)(^+ ^exp(2i((/)2 - 0i))) exp(2^C) C +00 
P\/7\/crexp(i0i + ^/c^Cl) C -00 

(4.10) 

^ / i^P^ exp(2^^i) + exp(2(^ + ^)Ci)) exp(2^C) C ^ +00 ,^ 
~ 1 1 + cxp(2ysr6) C ^ -00 



and hence we obtain the expansion 

(v^/pVT^) exp(i0i + ^{i)(y^+ y^exp(2i(02 - 0i))) 



C ^ +00 



l + (W7)exp(2^6) (4.12) 

( i + p2^exp(2ysrei) ^ °° 

where C — ^ ±00 corresponds to t ^ ±00. In the asymptotic past, this expansion (4.12) has 
the form of a 1-soliton solution 



~ ^ — ^ = lii , i ^ -00, = const. 4.13 

1 + exp (2^01(^1 -aj) 

in which the moving coordinate ^1 is shifted by 

ar = -ln(pV7)/^ (4.14) 

Similarly, in the asymptotic future, the expansion (4.12) again has the form of a 1-soliton 
solution 

t^+oc, 6= const. (4.i5) 

l + exp(2^(^i -aj)) 
in which the moving coordinate is now shifted by 

at = In(pv^)/^ (4.16) 
while in addition there is a phase shift given by 

^+^exp(2i(02-0i)) 

= exp(zi/i) (4-17) 

7 

where 

= arctan ( v^^;^2(0, - </>^) A (4.18) 
VVcT+Vc^ cos 2(01 -02)7 ^ ^ 
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Second we hold ^2 fixed and asymptotically expand F and G for large C = tAc — Ax with 
^1 = ^2 — C- After neglecting subdominant exponential terms, we obtain 

{-^^exp(i02 + A/C26) (\/c2 + ^/clexp{2i{(f)l - 02))) C ^ +00 
Vl V VV V ^^ .^g^ 

P^/l^/c2exp{i(f)2 + {^2 + 2^)C2) exp(-2^C) C ^- -00 



1 + 7p^ exp(2yc^^2) C +00 

(7p2 exp(2v^5) + exp(2(v^ + v^)^)) exp(-2^C) C -00 



(4.20) 



which yields the expansion 

\/^lpl\n) exp(#2 + v/c^6)(v^ + a/cT exp(2i(0i - ^2))) 



+00 



u 



l + p2^exp(2^{2) 

C ^ -00 



l + (l/p2^)exp(2V^6) 



where C, ±00 corresponds to t ^ ±00. In the asymptotic past, the expansion (4.21) has 
the form of a 1-sohton solution 



^exp 202 + 6 - 02 _ , , 

u ~ ^ — — = 1^2 , i ^ -00, ii = const. (4.22) 

1 + exp (2^(42 -^2)) 

in which the moving coordinate ^2 is shifted by 

= ln(pV7)/A/^ (4.23) 
In the asymptotic future, this expansion (4.21) similarly has the form of a 1-soliton solution 

l + exp(2^(^2 -aj)) 
in which the moving coordinate ^2 is now shifted by 

a+ = -ln(pV7)/Vc^ (4.25) 
while in addition there is a phase shift given by 

^+ycrexp(2i(0i-02)) . . . 
= exp(zi/2) (4.26) 

7 

where 

V = arctan { ^(02 - 0i) \ ^7) 

Va/^+a/^cos2(02-0i); ^ ■ ' 



Thus, for t — )■ ±00, the 2-soliton solution (4.4)-(4.6) asymptotically has the form of a 
superposition ~ itf + itf of a fast soliton itf and a slow soliton rtf , with speeds C\ and C2. 
As a consequence, this solution has the same conserved momentum (2.35), energy (2.36), 
Galilean energy (2.38), and center of momentum (2.39) as the mKdV 2-soliton solution. 
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4.2. Asymptotic position and phase shifts. 

In the 2-soliton solution u of the Sasa-Satsuma equation, the positions of the fast soliton 
uf and the slow soliton in the asymptotic past (t — )■ — oo) and future {t +00) are 
determined by the shifted moving coordinates = ~ = — cif- We thus see 

that the fast soliton ui retains its shape and speed, but gets shifted forward in position by 

. + - 1 , (v^ + V^)\/ci + C2 + 2^/5Ic^ COS 2(01 

A., = a,-a,=^ln^ (V^-^)^ )>' ^'-''^ 

while the slow soliton U2 similarly retains its shape and speed, but gets shifted backward in 
position by 



. + - 1 1 ( \/cr + yc^) + C2 + 2^5ic^cos 2 - 02 , 

Ax2 = - 02 = —In — < 4.29 

\/c2 \ (v/cT-x/c^)^ 



These asymptotic shifts satisfy the relation 

^/^Axi + v^Ax2 = (4.30) 

which can be understood as a consequence of the motion of the center of momentum of the 
2-soliton solution u in the same way as for the Hirota equation. Interestingly, in contrast to 
coUisions of Hirota solitons, here the shifts (4.28) and (4.29) depend on the relative phase 
angle 0i — 02 between the fast and slow solitons in the collision. 

Even more interestingly, in the collision, both the fast and slow solitons undergo a shift in 
phase given by (4.18) and (4.27) respectively. The features of these asymptotic shifts can be 
understood from the reflection properties of the 2-soliton solution as follows. First we write 
this solution in the equivalent form 

N P\/7(v^exp(i0i)(exp(-6l2) + exp(ii/i + 62)) + v/c^exp(i02)(exp(6li) + exp(ii/2 - 6*1))) 

x] — — — ^ 

^' ' 2(cosh(ei + ^2)+p'7Cosh(ei-^2)+cos(0i-(/)2)((v/^+v^)p2-l)) 

(4.31) 

with 

01 = ^i{x-cit), ^2 = v^(a:-C2t), (4.32) 

where we have used a space-time translation x ^ x — Xq and t — )■ t — to to shift the centers 
of the moving coordinates (4.9) to the positions 

Xi= Xq — CiIq = 0, X2 = Xq — 02^0 = (4.33) 

The solution (4.31) of the Sasa-Satsuma equation then exhibits an invariance 

■u(— t, —x) = exp{ii')u{t,x) (4.34) 

where the phase factor is given by 

^exp(i0i) + ^exp(i02) , oa w (■( < oa w or\ 

ex.p{iu) = ^ ^ = exp(z(i/i + 20i)) = exp(z(i/2 + 202)) (4.35) 

7 

in terms of the asymptotic phase shifts ui and z/2. Next we shift the phases 

01 ^ 01 - 1^/2, 02 ^ 02 - 1^/2 (4.36) 
which corresponds to an overall phase rotation of the solution 

u ^ ex.p{—iiy/2)u — u (4.37) 
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This does not affect the relative phase angle 0i — 02 between the fast and slow solitons in 
the collision, whereby the resulting solution 

. . . pyFf{^cosh{02 + iiyi/2) + ^cosh(gi - 

''^'^ ~ (cOSh(ei + e2)+p'7COSh(^l-e2)+COs(0i-02)((y^+V^)p2-l)) ^" 

is invariant under space-time reflection x —x, t —t, combined with phase conjugation, 

u{-t,-x) = iL{t,x) (4.39) 

For t — >■ ±cxD, this 2-soliton solution u of the Sasa-Satsuma equation describes a collision 
where, in the asymptotic past (t — )■ —00), a fast soliton with speed Ci, phase = 
01 — i//2 = —nui/2, and center of momentum xr(^) = + ^i^ overtakes a slow soliton ■U2 
with speed C2, phase 0^ = 02 — = —1/2/2, and center of momentum X2 (^) = ^^2^ + ^2^- 
In the asymptotic future {t — )■ +00), the fast soliton Ui undergoes a shift in both position 
xti^) = + Cit and phase (pf = Ui + (t)~[ = z/i/2, while the slow soliton U2 similarly 
undergoes both a position shift xti'^) = ^2 + and a phase shift 02 = ^2 + 02 = ^2/2, 
where these phase shifts are related by 

0r-02- = -(0i^-0^) (4.40) 

due to the reflection property (4.39). 

This asymptotic phase relation (4.40) is equivalent to 

i/l - i/2 = -2(01 - 02) (4.41) 

Thus, surprisingly, the relative phase angle between the fast and slow solitons is not preserved 
in the collision but instead changes sign. 



4.3. Interaction profile. 

The invariance property (4.34) of the 2-soliton solution u shows that the separation in 
positions of the peak amplitude of the fast and slow solitons uf and uf in the collision will 
be a minimum at time t — when the amplitude |m| is an even function of x. The shape 
of \u\ Bit t — thereby deflnes the interaction profile of the 2-soliton solution u, which can 
be understood to be the moment of greatest nonlinear interaction between the fast and slow 
solitons in the collision. This profile is characterized by the convexity of |m(0, a;)| at a; = 0. 

By an explicit calculation, we find that the convexity is given by 

\u{0,x)\^^^^ = 

(x/cT - \/c^) V2ciC2(CciC2 + B,/c^{ci + C2) - Ajcj + cj)) 

^7(^(^ + ^ + 7) -L(^ + ^))(g(v^ + ^)- +V^-i)) 

where 

A = cos^ A0 + 3 cos A0 + 3 
B = cos^ A0 + 6 cos A0 
C = 18cosA0 + 10 
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(4.43) 
(4.44) 
(4.45) 



and 

K^ci + C2 + 2 cos((/)i - 02) (4.46) 

L = 2y^(l - cos'(0i - 02)) (4.47) 

P = ci + C2 + 4^0^ (4.48) 

Q = 2y^(4 + 3 cos(0i - 02)) (4.49) 
The sign of the convexity (4.42) depends only on the ratio of speeds and the relative phase 

r = ci/c2, A0 = 01 — 02 mod 27r (4.50) 

of the fast and slow solitons. Since we have Ci > C2 > 0, tt > 0i > — tt, tt > 02 > — tt, these 
parameters are restricted to the respective intervals 

r > 1, TT > |A0| > (4.51) 

To determine the conditions under which the convexity is positive or negative, we will sep- 
arately consider the signs of the factors in the numerator and denominator expressions. 
The sign of the numerator in the convexity (4.42) is given by 

a = sgn{Cr + B^/^{r + l)-A{r^ + l)) (4.52) 

This expression can be factorized 

Cr + B^/F{r + l)-A{r^ + l) = A{{{D + B)/2A)y/^-r-l){{{D-B)/2A)^/^+r + l) (4.53) 

by means of the identity 

(D - B)(D + B) ^ 8>l(cos^ A0 + 12 cos A0 + 8) (4.54) 

where 

D = {3 cos A0 + 4) Vcos2 A0 + 12 cos A0 + 12 (4.55) 

We note this identity (4.54) also implies the inequality 

D> B (4.56) 

as follows. By evaluating D — B ior Acf) — arccos(— 6 + 2\/7), which is the sole real root of 
the right-hand side of the identity, we find D — B = 8{3\/7 — 7) is positive. This implies 
sgn{D — B) > holds for all values of A0 in the interval (4.51), since D — B may change 
sign only at a value where the identity (4.54) vanishes. 
To continue, from the inequality (4.56), we see 

{{D-B)/2A)^- + r + l>0 (4.57) 

holds throughout the interval (4.51), and hence the factorization (4.52)-(4.53) yields 

a ^sgn{{{D - B)/2A)^/^ + r + l) (4.58) 

which determines the sign of the numerator. 

The sign of the denominator in the convexity (4.42) is determined by 

S = sgn(2£;Vf ( vr + 1) - (r + 1 + A^r) ( + 1 - 7)) (4.59) 

where 

£; = 3cosA0 + 4 (4.60) 
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This sign can be evaluated in terms of the roots of the right-hand side of the identity 

{2E(r + V^) - (r + l + 4v/^)(v^+l -7))(2£;(r + v^) - (r + 1 + 4v/f)(v^ + 1 + 7)) 

= v^(4r(v^ - 1)^(1 - cos A0)^)^(A(r2 + 1) - (r + 1) - Cr) (4.61) 

There are two real roots contained in the interval (4.51), which are the only values of A0 
and r where 6 may change sign. For the root A0 = 0, we find d = sgn(14A/r(A/r + 1)) is 
positive since r > 1. Similarly, for the other root given by r + 1 + {{D — B)/2A)y/r — 0, we 
find S is again positive. Hence this implies 

5>0 (4.62) 

holds throughout the interval (4.51). We remark that the same argument can be used to 
show that the expression inside the square-root factor in the denominator of the convexity 
(4.42) is positive for all values of A0 and r in this interval. 

Therefore, the signs of the numerator (4.58) and the denominator (4.62) yield the sign of 
the convexity 

sgn(|M(0,x)|^^|^^p) = a (4.63) 
From expression (4.58) this sign is a quadratic polynomial in ^/r with two roots 

,/FZ={{B + D)/4A) ± ^/{{B + D)/4A)^-1 (4.64) 

which satisfy r+r_ = 1. The two roots are real and positive iff {{B + D)/AAf > 1. This 
condition can be expressed as > {4 A — B^ which simplifies to 

(6 cos A0 + 1) (cos^ A0 + 3 cos A0 + 3) > (4.65) 

Hence, the condition for r+ and r_ to be real and positive is cosA0 > —1/6. In this case 
the roots have the property r_)_ > 1 > r_ such that r_|_ = r_ = 1 when A0 = arccos(— 1/6). 
As a result, since r > 1, the convexity sign (4.63) is given by 

sgn(y^- ^/r) if |A0| < arccos(-l/6) ^ 0.557r 
< if \A<f)\ > arccos(-l/6) ^ 0.557r ^ ' ' 

where 

r+ = ^ + {B + VBiB + D)/2 + A{C-2A)) (4.67) 

Thus when the relative phase angle |A0| is greater than arccos(— 1/6) ~ 0.557r, the 2-soliton 
interaction profile |m(0, x) \ will always have a single peak at a; = 0, whereas when the relative 
phase angle jA^j is less than arccos(— 1/6) ?s 0.557r, the profile |ii(0,x)| will have either a 
single peak at a; = if the speed ratio is Ci/c2 > r+ or a double peak around a; = if the 
speed ratio is C1/C2 < r+. Additionally, for a double peak, the profile will always have an 
exponentially diminishing tail. For a single peak, the profile instead will have either a pair 
of side peaks around x — or just an exponentially diminishing tail, depending on whether 
I A0| is greater than or less than a certain critical value for which the profile exhibits a saddle 
point at some ,t 7^ as determined by the conditions |M(0,a;)|^^ = |-u(0,a;)|2. = (which we 
can solve numerically). The critical angle |A0| and the critical speed ratio Ci/c2 = r+ are 
shown in Fig. 13. 

In the case of a double peak profile, the interaction is a bounce-exchange, i.e. where the 
fast and slow solitons first bounce and then exchange shapes and speeds at t — x — 0. See 
Fig. 14. 
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Figure 13. Types of 2-soliton interactions for the Sasa-Satsuma equation 

In the case of a single peak profile without side peaks, the interaction of the fast and slow 
solitons is a merge-split, i.e. where the solitons first merge together at t = x = and then 
split apart. See Fig. 15 and Fig. 16. 

The case with side peaks is an absorb-emit interaction, i.e. where the slow soliton gradually 
is first absorbed by the front side of the fast soliton and is then emitted from the back side 
of the fast soliton. See Fig. 17 and Fig. 18. 

The special case of a saddle profile is shown in Fig. 19. 
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Figure 14. Sasa-Satsuma equation 2-soliton interaction with Ci/c2 = 3, 0i — 
02 = O.SStt 



M Iu| |u| 




Figure 15. Sasa-Satsuma equation 2-soliton interaction with Ci/c2 = 3, 0i — 
02 = O.SStt 



N |u| |u| 




Figure 16. Sasa-Satsuma equation 2-sohton interaction with Ci/c2 = 50, 
01 — 02 = 0.357r 
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Figure 17. Sasa-Satsuma equation 2-soliton interaction with Ci/c2 = 3, 0i — 
02 = 0.8557r 
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Figure 18. Sasa-Satsuma equation 2-soliton interaction with Ci/c2 = 50 
01 - 02 = O.SSStt 




Figure 19. Sasa-Satsuma equation 2-sohton interaction with Ci/c2 = 3, 0i — 
02 = 0.850327r 
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5. Comparison of soliton collisions for the Hirota and Sasa-Satsuma 

equations 



Up to a constant phase factor, the 2-sohton solution for both the Hirota and Sasa-Satsuma 
equations reduces to the mKdV 2-sohton solution when the relative phase angle A(/) between 
the fast and slow solitons in the collision is (same-orientation case) or ivr (opposite- 
orientation case). 

In the Hirota 2-soliton solution, the shift in positions of the fast and slow solitons is exactly 
the same as in the mKdV 2-soliton solution, which depends only on the respective speeds 
Ci and C2 of the solitons. Moreover, in both the Hirota 2-soliton solution and the mKdV 
2-soliton solution, there is no shift in the phase angles or orientations of the solitons. 

In contrast, the position shifts of the fast and slow solitons in the Sasa-Satsuma 2-soliton 
solution depend on their relative phase angle A0 in addition to their speeds Ci and C2. See 
Fig. 20. More remarkable is that the collision changes the sign of the relative phase angle. 

|u| 







' l-SoHtoa Hirota 2-Soliton Sasa-Satsuma 2-Soliton 



Figure 20. Sasa-Satsuma and Hirota 2-soliton interactions with Ci/c2 
<^>^ — (t>2 = O.Stt 



Both the Sasa-Satsuma and Hirota 2-soliton solutions describe three different types of 
collisions, which are separated by a critical speed ratio and a critical phase angle shown 
in Fig. 13 and Fig. 12. Interestingly, for collisions in the Sasa-Satsuma case, the angle 
|A0| ?a O.SSvr at which the critical speed ratio r_|_ approaches the limit ratio of 1 is strictly 
less than the critical phase angle |A0| ^ 0.712-7? when Ci/c2 approaches the same limit ratio 
of 1, whereas these angles are both equal to vr for collisions in the Hirota case. 

6. Concluding remarks 

Our work in this paper studying the interaction properties of complex mKdV solitons can 
be extended in at least three interesting directions. 

First, the Hirota equation (3.1) and the Sasa-Satsuma equation (4.1) are gauge-equivalent 
to a third order NLS equation [15] 



+ IqI'^q) + o:\q\^qx + (3\q\\q\xq + qx 
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(6.1) 



under the Galilean-phase transformation 

u{t, x) = q{i, x) exp ^±^.^^(5; — vi)^ (6.2) 

where v > is a speed parameter, with a = 24, /3 = in the Hirota case and a = /9 = 12 in 
the Sasa-Satsuma case. Our main results (cf. Sec. 5) concerning the properties of collisions 
of solitary waves for the Hirota and Sasa-Satsuma equations will thus directly carry over to 
solitary waves of the form 

/c 

q{i, x) — — exp (icj) + iut) sech {\/c{x — ci)) (6.3) 

with 



u ^±^/vJ2{c + v/2), c = c + 3^;/2>0 (6.4) 

for the corresponding cases of the NLS equation (6.1). In particular, such sohtary waves will 
exhibit three distinct types of collisions, which are separated by a critical speed ratio and a 
critical phase angle in terms of the speeds Ci, C2, the phases 0i, 02? and the frequencies Wi, 
U2 of the solitary waves in the collision. For all collisions, the waves undergo position shifts 
in both the Hirota and Sasa-Satsuma cases, as well as phase shifts in the Sasa-Satsuma case, 
such that the shifts depend only on the parameters ci, C2, and 0i — 02- 

Second, as indicated by the correspondence (6.2), the Hirota and Sasa-Satsuma equations 
each admit more general solitary wave solutions given by the form 

ti(t, x) = exp{iip{t))f{x — ct) (6.5) 

where 

i){t) = ujt + (t) (6.6) 

is a time- varying phase, with constant parameters w, 0, and where c > is the wave speed. 
Collisions of such waves should display interesting interaction properties with new features 
beyond those with a; = studied in Sec. 3 and Sec. 4. 

Second, both the Hirota and Sasa-Satsuma equations have a natural multi-component 
generalization given by the two known types of 0(iV)-invariant intcgrable mKdV equations 
[17] for a iV-component vector variable u. In particular, under the identification between a 
2-component vector u — (mi, U2) and a complex scalar u — Ui + iuii the vector version of the 
Hirota equation (3.1) is given by 

Ut + 2A\u\^u^^u^^^^Q (6.7) 
and the vector version of the Sasa-Satsuma equation (4.1) is given by 

Ut + 12|-u|^-Ua; -|- 12(-U ■ Ux)u + Uxxx = (6-8) 

For all > 2, these two vector equations are integrable and admit vector solitary wave 
solutions of the form 

u{t,x) = f{x -ct)ii) (6.9) 

with wave speed c > 0, where ip is an arbitrary constant unit vector and / is the sech 
solitary wave profile for the real scalar mKdV equation (2.1). In forthcoming work, we plan 
to generalize the results in Sec. 3 and Sec. 4 to study the interaction properties of collisions 
of vector solitons of both equations (6.7) and (6.8) for N > 2, where the collision involves 
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a fast soliton with vector orientation ■0i and speed ci overtaking a slow soliton with vector 
orientation ■02 and speed C2. 
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Appendix A 

The 2-soliton solution of the Sasa-Satsuma equation (4.1) can be derived most easily by a 
computational-ansatz version of the Hirota method [16] as follows. We first use the standard 
rational transformation 

u = G/F (A.l) 

with F real and G complex. This transformation converts the Sasa-Satsuma equation into 
the equivalent rational form 

= + Dl{G,F)) + -^^{GGDMG)) + j-,{2m' ' 3Dl{F,F))D,{G,F)) 

(A.2) 

as written in terms of Hirota's bilinear D operator given by 

D{a, b) = bDa - aDb (A.3) 

D^{a, b) = bD^a + aD% - 2{Da)Db (A.4) 

D^{a, b) = bD^a - 3{Db)D^a + 3{Da)D% - aD^ (A.5) 

Through the introduction of an auxiliary variable H, there is a natural splitting of equation 
(A.2) into a bilinear system of equations 

Dl{F, F) = 8GG, Dt{G, F) + Dl{G, F) = iHG, iFH = 6L>^(G, G) (A.6) 

for the real variables F, H and the complex variable G. 

We now observe that the 1-sohton solution of the bihnear system (A.6) is given by the 
simple exponential polynomials 

G = Aexp(e), F = 1 + Sexp(2^), H^Q (A.7) 

with 

e = kx-ut, u = k^, (A.8) 

where A and B are respectively complex and real parameters satisfying the algebraic relation 

= k^B (A.9) 

In particular, if we write A — \ A\ exp(i0) by a polar decomposition and express B — {\A\/kY, 
this yields 

G = kexp{k^), F = 1 + exp{2k^) (A. 10) 

in terms of the moving coordinate 

^ = x-kH-Xo (A.ll) 
which is centered at the initial position 

xo^Hk/\A\)/k (A.12) 
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Then the rational exponential solution obtained from expressions (A.l) and (A. 10) matches 
the solitary wave solution (4.2)-(4.3) having speed c — > and center of momentum 
X{t) = kH + xo. 

To obtain the 2-soliton solution of the bilinear system (A. 6), we use the ansatz 

G exp(^i) + ^2 exp(^2) + Ci exp(^i + 2^2) + C2 exp(^2 + 2^^i) (A. 13) 

F =1 + Siexp(2^i) + B2exp(2^2) + /exp(^i + ^2) + Jexp(2^i + 2^2) (A.14) 

and 

i/ = ii'exp(^i + ^2) (A. 15) 

with 

9i — kix — ujit, 02 — k2X — uj2t, (A. 16) 

where Bi,B2, 1, J, K are real parameters, and Ai,A2,Ci, C2 are complex parameters. In this 
ansatz, the form (A.13)-(A.14) for the variables G, F is motivated by the rational exponential 
form of the mKdV 2-soliton (2.13), while the form (A. 15) for the auxiliary variable H comes 
from balancing the highest-power terms in the third equation in the bilinear system (A. 6). 
Substitution of G, F, H into the system (A. 6) leads to three equations, which are each a 
polynomial in the exponentials exp(^i) and exp(^2)- The coefficients of these polynomials 
directly yield an overdetermined bilinear system of algebraic equations that can be solved 
for the parameters 5i, i?2, 1, J, K, Ci, C2 in terms of the pair of complex parameters Ai,A2. 
This gives us the result (obtained by a Maple calculation) 

uji = kl, 0J2 = kl, 9i = ki{x - kit), 62 = k2{x - kji), (A.17) 
I /I 12 I /I 12 

^ {ki - k2){kiAiA2 - k2A2Ai)A2 ^ (fci - k2){kiAiA2 - k2A2Ai)Ai 

' ki{k, + k2y ' ' ki{k, + k2y ' ^ ^ ^ 

_ , A1A2 + A2A1 j_ {ki- k2y{kiA2Ai - k2AiA2){kiAiA2 - ^2^2^) . 
{k, + k2r ' kfkl(k, + k2r ' ^ ^ 

K ^ 6{ki - k2){iA2Ai - iAiA2) (A.21) 

Note this solution (A.13)-(A.20) contains a pair of arbitrary phases and positions corre- 
sponding to the parameters 

Ai = 1 exp(iV'i), A2 = 1^2 1 exp(iV'2) (A.22) 

We now simplify the form of the 2-soliton solution (A.13)-(A.14) and (A.17)-(A.20) by 
the following steps. First we observe 

klCi/A2 = klC2/Ai = l^ill^lC (A.23) 

J= (CiM2)(C'2/Ai) = {\A,\\A2\\C\y/{k,k2f (A.24) 

where 

^ = 7r^—rlikiexp{tA^) - fc2exp(-^A^)), = ^1 - ^^2 (A.25) 

(fci + k2y 

Next we introduce the moving coordinates 

$,1 — X — kft — xi, ^2 — X — k^t — X2 (A. 26) 
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centered at initial positions 

x^ = In {h/{\A,\^/\C\))/h, = In {k^/ {\A2\^/\C\)) 1^2 (A.27) 

which are determined in terms of |Ai| and \A2\. Then the exponential polynomials (A. 13) 
for G and (A. 13) for F are given by 

G = fA:iexp(gi)(exp(i'^i) + exp(i('j/;2 + ip)) exp(26'2)) 

A;2exp(^2)(exp(iV'2) +exp(i(V'i +V'))exp(2^i))) (A.28) 
F=l + exp(2^i + 2^2) + p (exp(2^i) + exp(2^2) 

+ ^kik2{ki + A;2)"^(exp(iAV') + exp(-iAV')) exp(ei + ^2)) (A.29) 

in terms of 

= A;i6, 92 = k2i2 (A.30) 

and 

exp(z2V') = (7/(7 = 1 7 — -r-T ; . (A. 31) 

^ ^ ' A;iexp(-2AV^) - A;2exp(?A^) ^ ' 

These expressions (A.28) and (A.29) match the form of the colliding solitary wave solution 
(4.4)-(4.6) if we write 

V' = 02-0i, = V'2 = </>i-AV' (A.32) 
and use relations (A. 31) and (A. 25) to express 

exp(iA'0) — — (/ci exp{iip) + k2exp{—iip)) — exp(— i(z/i — -0)) = exp{—i{i'2 + (A. 33) 

and 

1^1 = ,,^^7,%^^"^^' = —2 (A.34) 
(fci + k2y 7 

where 

7 = Jkl + kl + 2k,k2Cos{2ij), p = ^j5l±5 (A.35) 
V ki — k2 
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